Quantum States of Neutrons in Magnetic Thin Films 
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We have studied experimentally and theoretically the interaction of polarized neutrons with mag- 
netic thin films and magnetic multilayers. In particular, we have analyzed the behavior of the critical 
edges for total external reflection in both cases. For a single film we have observed experimentally 
and theoretically a simple behavior: the critical edges remain fixed and the intensity varies accord- 
ing to the angle between the polarization axis and the magnetization vector inside the film. For the 
multilayer case we find that the critical edges for spin up and spin down polarized neutrons move 
towards each other as a function of the angle between the magnetization vectors in adjacent ferro- 
magnetic films. Although the results for multilayers and single thick layers appear to be different, in 
fact the same spinor method explains both results. An interpretation of the critical edges behavior 
for the multilyers as a superposition of ferromagnetic and antifferomagnetic states is given. 

PACS numbers: 74.78.Fk, 61. 12. Ex, 75.25.+Z 
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INTRODUCTION 

Neutron reflectivity goes back to 1946 when it was 
first used by Fermi and Zinn for the determination of 
coherent scattering lengths 0. Subsequently inserting 
a polarizer and analyzer to produce a polarized neutron 
beam, Hughes and Burgy started to perform polarized 
neutron experiments as early as 1951 Q{. It was fore- 
seen by the authors of this work that beams of completely 
polarized neutrons will be useful in the study of mag- 
netic and nuclear properties. The next major achieve- 
ment of Polarized Neutron Reflectometry (PNR) was the 
prediction of spin-flip reflectivity by Ignatovich (1978) Q 
and the pioneering experiments on magnetic surfaces by 
Felcher (1981) While specular PNR is widely rec- 
ognized as a powerful tool for the investigation of mag- 
netization profiles in magnetic heterostructures the 
description of off-specular scattering from magnetic do- 
mains is still under development ;8j. In spite of these 
important developments there is still a confusion concern- 
ing the quantum states of neutrons in a magnetic sample. 
Here we show unambiguously that the neutron has to be 
treated as a spin 1/2 particle 0,0 m eacn homogeneous 
magnetic layer. This is at variance with the conventional 
description of neutron reflectivity, which often considers 
the neutron magnetic potential as a classical dot prod- 
uct 0, [13 El 03. 

Neutrons interact with a magnetic thin film via the 
Fermi nuclear potential and via the magnetic induction. 
Thus, the neutron - film interaction hamiltonian includes 
both contributions: V = V n + V m = (fr 2 /2m)AirNb- \iB, 
where m is the neutron mass, N is the particle density 
of the material, b is the coherent scattering length, |/x| 
is the magnetic moment of the neutron, and \B\ is the 
magnetic induction of the film. Unconventionally, how- 



ever, neutron reflectivity treats the dot product between 
the magnetic induction and neutron magnetic moment 
classically: V m = fiB = ±\fi\\B\ cos(0), where 9 is the 
angle between the incoming neutron polarization direc- 
tion and the direction of the magnetization inside the 
film. Writing the magnetic potential as a classical dot 
product implies that the neutron energies in the mag- 
netic layer have a continuous distribution from - to 
+ This predicts that the critical angle for total 

reflection depends on the angle between the direction of 
polarization and the direction of the magnetic field inside 
the layer: 
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where a is the glancing angle to the surface, A is the 
wavelength of the neutrons, and is the critical scat- 
tering vector. There are experimental data |lll 1 1 3| on 
magnetic multilayers which apparently confirm this be- 
havior. Therefore, the classical representation appears 
to provide a convenient and trans pare nt way to describe 
the experimental observations [9I lift Ull Hijl . 

From the Stern-Gerlach experiment we know that 
there are only two eigen states for the spin 1 /2 particles 
in a magnetic field. Therefore, the eigen wave number 
of a neutron in a magnetic thin film has two proper val- 
ues. After solving the Schrodinger equation one obtains 
two eigen wave numbers for neutrons in a magnetic film: 
k± = TpKKj ± |/j||B|). They correspond to two possi- 
ble states of spin orientation: one for the case, when the 
spin is parallel to the magnetic induction, and the other 
one for the antiparallel orientation. It follows that there 
are only two possible energies and consequently only two 
values for the index of refraction corresponding to the 
spin-up and spin-down states of the neutrons. Therefore, 
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QM predicts that there are only two critical angles for 
the total reflection: one corresponding to the R + and 
one to the R~ reflectivity 
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Obviously there is a contradiction between the quantum 
mechanical prediction (Eq[2J) and the prediction based 
on the classical representation of the magnetic poten- 
tial(Eq^|: quantum mechanics predicts that the spin 
states of the neutron is determined by the magnetic in- 
duction in the sample, whereas classical representation 
of the magnetic potential, supported by experiments on 
magnetic multilayers, assert that the spin states of the 
neutrons is fixed by the incident polarization axis. 

Here we describe an experiment which provides direct 
and unambiguous evidence for the spin states of neutrons 
in magnetic media. The goal is to find a system where 
the angle between the neutron polarization and direction 
of the magnetization inside of the film can be fixed and 
controlled. Then we measure the R + and R~ reflectiv- 
ities and determine whether the position of the critical 
edges changes as a function of the angle 6, or whether the 
critical edges stay fixed, and only intensity redistributes 
between reflections R + and R~ with change of 9. The 
easiest way to control the angle 9 is to rotate the mag- 
netic film and therefore the magnetization direction with 
respect to the neutron spin polarization, which remains 
fixed in space outside of the sample. This requires that 
the film should have a high remanent magnetization. Ad- 
ditionally, the film thickness should exceed the average 
neutron penetration depth |l5j . The last requirement 
is essential in order to avoid neutron tunnelling effects 
which will hinder the precise determination of the criti- 
cal edges. 



SAMPLE CHARACTERIZATION BY MOKE 
AND PNR 



To fulfill the aforementioned requirements, we have 
chosen a 100 nm thick polycrystalline Fe film deposited 
by rf-sputtering on a Si substrate. The thickness of the 
Fe films was about 4 times larger than the average pen- 
etration depth 1 / y/ '2mVjv / ■ The Fe film was covered 
with thin Co and CoO layers, the latter one protecting 
the Fe film from oxidation. For sample characterization 
at room temperature we first recorded hysteresis loops 
with the magneto-optical Kerr effect (MOKE). A series 
of hysteresis loops were taken with the field parallel to 
the film plane but with different azimuth angles of the 
sample. A typical hysteresis loop is shown in Fig. ^ 
The coercive field is about 20 Oe and the remanence is 
high. A plot of the ratio between the remanent magneti- 
zation and saturation magnetization M rem /M sat versus 
the rotation angle about the sample normal is shown in 
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FIG. 1: Bottom: hysteresis loop of the polycrystalline Fe/Si 
sample measured by MOKE. Top: the behavior of the rema- 
nent magnetization as a function of the rotation angle ex- 
tracted from hysteresis loops. 
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FIG. 2: (Color online). Top: Polarized neutron reflectivity 
curves R + (solid black circles) and R _ (open black circles) 
of the Fe/Si sample. The black line is the simulated R + re- 
flectivity and the red (light gray) line is the simulated R~ 
reflectivity. The applied magnetic field was 2000 Oe. In the 
inset, the magnetic profile obtained from fitting the data is 
shown. Bottom: Experimental (open black symbols) and sim- 
ulated (black line) spin asymmetry ((R + — R~)/(R + + R~)) 
are plotted for the same sample in saturation. All lines in the 
top and bottom panels are fits to the data points using the 
GMM (for more details see text). The abscissa is the wave 
vector transfer: Q — 47rsin(a)/A. 
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Fig. ^ We conclude that the system has no macroscopic 
anisotropy and the remanent magnetization is 97.5% of 
the saturation magnetization. 

Neutron reflectivity experiments were performed us- 
ing the angle dispersive neutron reflectometer ADAM in- 
stalled at the Institut Laue-Langevin, Grenoble, which 
operates at a fixed wavelength of 4.41 A. The R + and Rr 
reflectivities and the spin asymmetry (R + - R~)/(R + + 
R~) in a saturation field of 2000 Oe are plotted in the 
top and bottom panel, respectively, of Fig. |2 The solid 
lines are fits to the data points using the PolarFit code 
based on the general matrix method (GMM) 7]. The 
fit and sample parameters are listed in Table ^ In or- 
der to obtain a high confidence of the fit parameters, all 
reflectivities were fittted together and with the same pa- 
rameter set. In general it is useful to fit first the spin 
asymmetry, for which geometrical and normalization pa- 
rameters drop out. 



TABLE I: Parameters of the Fe/Si sample obtained by fitting 
to the R + an R~ data shown in FigEl d is the layer thickness, 
a is the rms roughness, SLD is the scattering length density, 
and B is the magnetic induction in the ferromagnetic films. 
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From Table [H and the magnetic characterization 
(Fig. ^) we conclude that the sample fulfills the require- 
ments as concerns thickness, anisotropy, and rcmancncc 
as required for our experiment. 




FIG. 3: (Color online). Experimental (symbols) and simulated (lines) reflectivity curves (R+ and R") from Fe(1000 A)/Si 
sample. The abscissa is the wave-vector transfer. The two sets of R + (solid black symbols) and R~ (open black symbols) 
reflectivity curves were measured for four different 6 angles between the neutron polarization along Bo and the direction of 
the magnetic induction (B) which lies in the sample plane. The guiding field is Bo ~ 10 Oe. The blue (thick dark gray) lines 
are the simulated R+ reflectivities and the red (thin light gray) lines are the simulated R- reflectivities. In the right side the 
experimental geometry is shown. The figure shows that the critical edges Q\ and Q'L are not sensitive to the 9 angle. 



ROTATION EXPERIMENT 

The rotation experiment was performed as follows: the 
Fe layer was magnetized parallel to the neutron polariza- 
tion direction and then the magnet was removed. A small 
guiding field (H c ) is still present at the sample position 
in order to maintain the neutron polarization. Subse- 
quently a series of R + and R~ reflectivities, shown in 
Fig. 13 were measured for several in-plane rotation an- 
gles of the sample. We observe two characteristics of the 



reflectivities: (1) the critical edges arc fixed and indepen- 
dent of the in-plane rotation angle 9 (see Eq. |2J), and (2) 
the R~ intensity continuously increases at the expense 
of the R + intensity as a function of the 9 angle. Using 
the parameters obtained from the fit to the saturation 
data (see Table HJ and using the rotation ang le 9 as set 
during the experiment, we have simulated 14] the reflec- 
tivities using the GMM approach 7] which transparently 
predicts the behavior of the critical edges described by 
Eq.[3 The simulated curves are plotted together with the 



4 



experimental data in Fig. [3] There are no free parameters 
for these simulations, providing an excellent description 
of the experimental results. The fixed critical edges Q c + 
and Q c _ can easily be interpreted in the context of the 
neutron spin states in homogenous magnetic media as 
discussed in the introduction. 



EXPERIMENTAL DETERMINATION OF 
MAGNETIZATION ORIENTATION 

The sensitivity to the in-plane rotation angle of the 
magnetization is seen very clearly seen in the reflected 
intensities R + and R~ plotted in Fig. . It has been 
shown theoretically 0] that, for a single magnetic layer, 
the normalized spin asymmetry (nSA{9)) is directly re- 
lated to the 8 angle through the following expression: 



" s "< e > - lm =cos<9) 



(3) 



Now, we use our experimental data shown in Fig. [3] to 
confirm the validity of this equation. In Fig. 0] is shown 
the experimental normalized spin asymmetry and the co- 
sine of the experimental angles. The agreement between 
the experimental normalized spin asymmetry (symbols) 
and the cosine of the 6 angles (lines) set during the ex- 
periment is excellent over the whole wave vector transfer 
range. It follows that the magnetization orientation of 
a single magnetic layer with respect to the neutron po- 
larization outside the layer can be easily extracted ex- 
perimentally using Eq. For more complicated sys- 
tems a numerical fitting is still necessary. The nSA is 
an important measure of hysteresis loops. It was shown 
in Ref. [16j that nSA can be written, generally, as: 
nSA = Mu/M sa t for both, magnetization reversal via 
coherent rotation and via domain wall movement. This 
implies that nSA reproduces the hysteresis loops as mea- 
sured by SQUID or MOKE. Here we confirm the validity 
of the nSA for determining the magnetization reversal 
via coherent rotation. We mention that this equation is 
valid for samples which contain a single magnetic layer. 
Comparing MOKE or SQUID hysteresis loops with nSA 
is a very useful tool for the evaluation of magnetic do- 
main state and/or a reduced magnetization within the 
layer. 

MULTILAYERS 

Our next topic is to investigate the neutron spin states 
in multilayers with noncollinear magnetization of adja- 
cent layers. We simulated the reflectivity profile of a 
[Fe(60A)/Cr(8A))]4 /Si superlattice, with thicknesses of 
the Fe and Cr layers which are typical for many real 
superlattices 01 . For the simulation we used the free- 
ware code PolarSim 0] based on GMM for the calcula- 
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FIG. 4: (Color online). Solid black symbols: The experimen- 
tal normalized spin asymmetries (nSA(0) = S A(9) / (S A(0))) , 
plotted as a function of the wave- vector transfer Q. Lines: The 
three lines are the cosines of the corresponding angles set dur- 
ing the experiment. From top: cosine of 48° (thin gray (red) 
line) , 68° (thick gray (blue) line), and 88° (thin black line), 
respectively. The angles are the experimental rotation angles 
6 used during the experiment shown in Fig. |H] The experi- 
mental normalized spin asymmetries are assembled from the 
R + and R~ reflectivities shown in Fig [3] This figure shows 
that the equation nS A(9)=cos(9) is valid over the whole Q 
range for a single magnetic layer. It can be used to extract 
the angle 9 directly from the experimental reflectivities. 



tion of the reflection and transmission coefficient together 
with a full quantum mechanical description of the spin 
states [3 • In the simulation the choice of a Si substrate 
has the advantage that it does not obscure the critical 
edge of the (-) neutron state. In the top panel of Fig.[S]we 
show simulations of R + and R reflectivities for three an- 
gles 7 between the magnetization vectors of adjacent Fe 
films: 7 = (or ferromagnetic alignment); 7 = 100°; and 
7 = 170° (close to antiferromagnetic alignment). Our 
focus is on the behavior of the critical scattering vec- 
tor for total reflection. We observe that for 7 = the 
(+) and (-) critical scattering vectors are well separated 
and that they contain information about the saturation 
magnetization. When the 7 value increases, the critical 
edges approach each other. For an angle 7 = 180° (not 
shown here) there is no difference between the R+ and 
R- reflectivities. The main result from this simulation is 
the observation that the separation of the critical edges 
is a continuous function of the angle 7 between the in- 
plane adjacent magnetization vectors. The critical edge 
positions satisfy the following relation: 
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where is an effective nuclear potential. Clearly, for 

this geometry the angles 9 and 7/2 coincide if the neutron 
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polarization is parallel to the average field (B\ + Z?2)/2, 
where B\ and Bi are the magnetic field inductions of 
adjacent layers. Therefore, numerically, the eqs. f and 
3 are almost identical. However, there is a fundamental 
difference: similarly to the single layer, the 9 angle does 
not influence the position of critical edges, whereas the 
7 angle is solely responsible for the continuous shift. 

To shed more light on how the 9 and 7 angles affect 
the critical edges for polarized neutron reflectivity at the 
multilayers we simulated numerically the rotation exper- 
iment performed on the single layer. For a fixed coupling 
angle of 7 = 90, as it can be achieved also experimentally 
via biquadratic exchange coupling, the reflectivities R + 
and R~ are plotted as a function of the 9 angle. Here the 
9 angle is the angle between the incoming neutron po- 
larization and the direction of the average magnetization 
vector of two adjacent ferromagnetic layers. The results 
are shown in the bottom panel of Fig. [SJ We observe a 
similar behavior of the critical edges and intensities as for 
the single layer. While the positions of the critical scat- 
tering vectors Q c + and Q c _ remain fixed for a constant 
coupling angle 7, the intensity increases on the ex- 
pense of the R + intensity with increasing 9 angle. With 
this simulations we lift the contradiction stated in the 
introduction by showing that Eq. ^ is a particular case 
of Eq. 0] which, in turn, is in agreement with the QM 
description of the neutron spin states in magnetic media 
Eq. [21 The different behaviors of the critical edges for 
the case of a single homogeneous ferromagnetic layer and 
for a multilayer with alternating directions of the layer 
magnetization vectors now becomes obvious: in the mul- 
tilayer the neutrons are affected by an average magnetic 
potential which depends on the relative orientation of the 
magnetic induction in the individual layers. However, in 
both cases, single film as well as multilayer, the magnetic 
potential of the individual layers (V m — \fi\\B s \) enters 
the algorithm for calculating the reflectivities. 

It should be noted that the dependence of Q c + and Q°_ 
on the angle 7 in a multilayer is a general property of 
the periodic potential with different field orientation and 
magnitude. It is natural to expect that such a sample is 
a noncollincar fcrrimagnet with ferromagnetic field B f = 
(£?!+£? 2 )/2, 

B f = \B 1 + B 2 \/2 = Bcos(~//2), (5) 

and with antiferromagnetic field B a f = (B\ — £?2)/2, 

B af = \Bi - B 2 \/2 = Bsin(7/2). (6) 

Then, the critical edges can be expected to be given by 
Eq.H 

To further stress the origin of the effective nuclear po- 
tential V% term in Eq. |2J, let us consider the critical 
edge for non-polarized neutrons when scattered at a [Fe(x 
A)/Cr(y A))]oo multilayer. Naively, we may expect that 
the critical edge to be given by the Fermi interaction po- 
tential of Fe as it is higher than the potential of Cr. This 
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FIG. 5: (Color online) . Top: Simulation of polarized neutron 
reflectivities [R + (dashed lines, shifting from right to left) , R~ 
(solid lines, shifting from left to right)] for a Fe/Cr multilayer 
as a function of coupling angle (7) between the magnetization 
vectors of adjacent Fe layers. Bottom: Simulations of R + 
(solid lines, shifting from top down) and R~ (dashed lines, 
shifting from bottom up) as a function of rotation angle 9 for 
7 = 90° 



is, however, not the case. For a finite thickness x of the 
Fe layer and zero thick Cr layer, indeed the critical edge 
is equal to the critical edge of a single thick Fe layer. 
Vice versa, for zero thickness of Fe layer and finite thick- 
ness y for the Cr layer the critical edge is given by the 
Fermi potential of Cr. However, when both layers have 
finite thicknesses the critical edge of the multilayer will 
vary from the value for pure Fe to the value for pure Cr. 
Therefore, the critical edge for non-polarized neutrons re- 
flected from a multilayer not only depends on the Fermi 
potential of the two separate layers, but also on their 
individual thicknesses. 
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CONCLUSIONS 

In summary, we have analyzed the behavior of the criti- 
cal scattering vectors Q c + and Q c _ for total external reflec- 
tion of a polarized neutron beam for the case of homoge- 
neous ferromagnetic films and for antiferromagnetically 
coupled multilayers. For a single film we have observed 
experimentally and shown theoretically that the critical 
edges do not change as a function of the angle between 
the neutron polarization and the direction of the mag- 
netic spins inside the film. They fulfill the relation Eq.|2 

Qf = (V n ±\fi\\B a \), which directly reflects the 

spin states of the neutron beam in magnetic thin films. 
For multilayers we found that the critical edges for total 
external reflection move towards each other as a function 
of the coupling angle. Their position is well reproduced 

by the Eq. m Qt = J%£(V° ff ±\n\\B a \cos(>y/2)). 
The cos(7/2) dependence is not related to the neutron 
spin states in the magnetic media, but it is the result of 
the presence of a ferromagnetic field direction along the 
average field in the noncollinear ferrimagnetic. By choos- 
ing a fixed coupling angle 7 between the magnetization 
vectors of adjacent layers and rotating the sample, the 
critical edges behave again in accordance with the neu- 
tron spin states in homogeneous magnetic media. Practi- 
cally, the coupling angle in non-collinear superlattices can 
be inferred directly from the experimental data through 
the separation of the critical edges. For a single layer the 
orientation of the magnetization can be extracted exper- 
imentally from the spin asymmetry. 
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